Quantum phonons and the charge density wave transition temperature: a dynamical 

mean field study 
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We use the dynamical mean-field method to calculate the charge density wave transition temper- 
ature of the half-filled Holstein model as function of typical phonon frequency in the physically 
relevant adiabatic limit of phonon frequency Q much less than electron bandwidth t. Our work is 
the first systematic expansion of the charge density wave problem in Q/t. Quantum phonon effects 
are found to suppress T co severely, in agreement with previous work on one dimensional models and 
numerical studies of the dynamical mean field model in the extreme antiadiabatic limit (fi ~ t). We 
suggest that this is why there are very few CDW systems with mean-field transition temperatures 
much less than a typical phonon frequency. 
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I. INTRODUCTION 

The present-day understanding of electron-lattice in- 
teractions injnetals is based on the Migdal-Eliashberg 
(ME) theorytl, which exploits the smallness of typical 
phonon frequencies, f2, and temperatures, T, relative to 
typical electronic energies, t. This theory may be re- 
garded as an expansion in powers of the adiabatic pa- 
rameter 7 = max(f2,T)/t, in which only the leading non- 
trivial term is retained. In many cases the leading term 
is of order 7 and thus within ME theory most quantities 
of physical interest are independent of ion mass. A no- 
table exception is the superconducting transition temper- 
ature T c , which is controlled by a logarithmic divergence 
with an upper cutoff set by the phonon frequency, and 
involves an interaction (the Coulomb pseudopotential) 
which has a logarithmic suppression with a lower cutoff 
set by the phonon frequency. The resulting, rather com- 
plicated, dependence of T c on ion mass was determined 
by McMillarJi 

A series of recent experimental papers reporting iso- 
tope effects (i.e. ion-Haass dependence) for a variety of 
electronic propertiestiu motivates re-examination of this 
issue. Of particular interest here are reports of a pro- 
nounced increase of the charge ordering transition tem- 
peratures of various members of the 'CMR manganitc' 
family oLmanganese perovskites when 18 is substituted 
for 16 0. Era This leads us to re-examine the theory of 
charge density wave (CDW) instabilities in order to de- 
termine how the transition temperature T co is affected 
by ion mass. A theory applicable to materials of current 
interest must go beyond the Migdal-Eliashberg approxi- 
mation, must apply to the case of d — 3 spatial dimen- 
sions and must be able to treat the physically relevant 
limit 7 -C 1. 

The ^-dependence of T co has been studied. In one 



spatial dimension pioneering numerical work of Hirsch 
and FradkinE2l followed by analytical studies by Bour- 
bonnais and CarorO found a strong suppression of the 
zero temperatj«e-CDW order parameter with increasing 
f2. In d — 2t3i3 and d — 00EJ the quantum phonon 
problem has been attacked by a direct numerical ap- 
proach using quantum Monte Carlo (QMC) methods. 
However, numerical issues related to the mismatch be- 
tween phonon and electron frequency scales however re- 
stricted the QMC studiealdilJ to the antiadiabatic limit, 
Q ~ t. Also, in Refs. [T^-p^ the isotope effect on the 
CDW transition temperature was not studied in detail, 
although values for T co for a given phonon frequency 
were calculated. In d = 00 and for weak coupling 
strengths the CDW transition temperature is strongly re- 
duced with decreasing ion mass and a reduction ratia.of 
T co {Q = OM)/T CO (Q = 0) w 0.25 has been reportedMlia 

Especially in d = 00 considerable effort has been ex- 
pended to construct weak-coupling perturbation theories 
which continuously connect the adiabatic to the antia- 
diabatic limit oii3 One important outcome of this work 
was the suggestion that future research should focus on 
generalizations of the ME theory that work with dressed 
phonons (renormalized ME) but include higher order 
non-adiabatic effects such as vertex corrections. Our 
method is such an approach.— A large number of stud- 
ies based on the Lang-Firsovlla method have been pub- 
lished (see e.g. Ref. |T^). For one dimensional systems 
the Lang-Firsov results are in agreement with previous 
work.EH However, this is an uncontrolled approximation 
optimized for the single polaron problem and its applica- 
bility to metallic densities is not established. 

In this paper we construct an essentially analytic the- 
ory of the charge density wave transition temperature in 
the physically relevant limit 7^1. We build on previ- 
ous work which showed how to construct an expansion in 
7 within the dynamical mean-field approximations and 
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which analysed the classical limit (CI — 0) of the charge 
density wave problem.B Here we show how to calculate 
T co including the first non-Migdal term in 7 for arbi- 
trary values of Q/T co . We find that T co depends very 
strongly on CI, becoming very small as CI is increased 
(at fixed coupling) beyond the O = value of T co , but 
that the specific behavior depends strongly on whether 
spin-1/2 or spinless electrons are considered. For spinless 
electrons, relevant to half-metals such as CMR mangan- 
ites, a divergent isotope exponent is obtained in the weak 
electron-phonon coupling limit. 



II. MODEL AND FORMALISM 

In this paper we study the simplest model of electrons 
interacting with phonons, namely the Holstein modelEil 
-HW = H Q \ + H p h + i/oi-ph with 

ij i 

(3) 
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where we have absorbed the electron-phonon coupling 
into the phonon coordinate r which thus has dimension 
of energy as does the phonon stiffness parameter A. We 
further specialize to a mean electron density per spin 
direction of n = 1/2; this implies p, = 0. We also assume 
a bipartite lattice, which for our purposes we define as 
a lattice possessing a dispersion (Fourier transform of 

U-j) and a wavevector Q such that eg_g 



eg for all k 



in the Brillouin zone. An equivalent definition is that the 
lattice may be divided into two sublattices A and B, such 
that ti-j connects A sites only to B sites. The important 
feature of the dispersion is the bare density of states per 
site per spin p(e) — -h ^2j:5(e— eA. In this paper we use 
the semicircular form p(e) = -^^y/At 2 — e 2 . It is useful 
to introduce the dimensionless electron-phonon coupling 
A and adiabatic parameter 7 by 



(4) 



Here po = 1/irt denotes the bare density of states at the 
chemical potential. 

We now establish notation by reviewing the deriva- 
tion of the T co equation within dynamical mean-field the- 
ory (DMFT). The fundamental assumption of DMFT is 
the momentum independence of the electron self energy 
£„ = S(iw„) (u> n = 2-kT (ti+1/2)). Consequently, in the 
absence of spatial symmetry breaking, £„ can be derived 
from an effective impurity model specified by the action 

S({r nm }; {c„}) = (T/2A) £ r nQ (u 2 n0 /Cl 2 + l) r Qn 



r n0 - (2s + 1) Tr In [c„ S nm - T r,,, 



(5) 



In this action r nm are bosonic fields describing the 
phonons and c„ mean field functions describing the 
fermions (which have been integrated out). We have in- 
dexed the bosonic fields by the difference of two fermionic 
Matsubara frequencies uj nm = uj n — uj m . The factor 
(2s + 1) takes the values one and two for spinless (s = 0) 
and spin-1/2 (s = 1/2) electrons, respectively. The par- 
tition function is the functional integral over the bosonic 
fields 



Z = 



J V[r] exp(-S) 



(6) 



and is a functional of the effective field c alone. It is 
useful to define the impurity Green's function Q and self 
energy S by 



Gn = 



1 5 In Z 



1 



1 



2s + 1 S c„ 



2s + 1 Cn - S„({c„}) 



(7) 



The effective field is fixed by equating the local Green's 
functions of the original lattice model and the effective 
impurity model. For a semicircular density of states one 
obtains 



On = iuj n + p - t 2 Q n ({c n }) 



(8) 



For temperatures below the CDW transition the self 
energy acquires two components for the two sublattices 
A and B which calls for two effective impurity models 
given by the effective fields a and b respectively. The 
resulting self-consistency equations are: 



a n = iu) n + p-t 2 Q n ({&„}) 

b n = iu) n + p-t 2 Q n ({(!„}) 



(9) 
(10) 



Near the transition (assumed to be second order) we can 
expand a n — c n + e' n and b n — c n — e' n leading to e' n 



Jm -q^- c m . Making use of Eq. (0) we obtain 
1 



e„ = 



1 + (c n - £„) 2 A 2 ^ do 



(11) 



It is amvenient to rewrite Eq. (^Tj) in the more symmetric 
forrrJj'Ea 



v 1/2 e 

ran Am c r 



(12) 



Here we have introduced \ n = —Q n j(iu: n — £„), the ir- 
reducible two particle vertex T nm = J2 a > ®^n/®@m an( l 
e n = t' n Xr/ 2 ■ Eq. dl^) is the statement that the ma- 



trix -Xn^liiinXm^ ^as eigenvalue one and T co is the 
temperature at which this eigerpalue problem has a so- 
lution. It can be easily shownlij that \n is the 'bare' 
CDW susceptibility of the original lattice model x„ = 
-1/N J2k G n(k + Q) G n (k). 
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III. SOLUTION IN THE ADIABATIC LIMIT 

As written, Eq. ( |l2| ) is general, but not very useful 
because not only must the matrix be diagonalized nu- 
merically, but also determination of the matrix elements 
for arbitrary A and 7 requires performing the functional 
integral over phonon fields numerically (i.e. via quantum 
Monte Carlo) and this becomes prohibitively lengthly for 
physically relevant regimes of f2 and T. Here we turn 
this difficulty to advantage, exploiting the smallness of 
n and T to construct an analytic expansion in powers 
of 7. The expansion builds on the observations that the 
biggest contribution to a diagram (e.g. for the electronic 
self energy or the vertex) with n phonon loops is of order 
[max(T, fl)/t] n . This can be seen as follows. A generic 
diagram with n phonon loops contains n inter- 

nal sums over bosonic Matsubara frequencies u>i, . . . , u> n 
(here to n = w„o) each associated with one (or more) free 
phonon propagators given by 



(a) 



T>° = -A 



o 2 



n 2 



(13) 



If F(iu>i, . . . , iu> n ) depicts the appropriate product of 
Green's functions and phonon propagators we can write 
D (n) = T n J2i n F(iuj 1 ,...,iuj n ). In the classical 
limit, Q/T — > 0, the free phonon propagators become 
delta functions and therefore D^ n ' ~ T n as noted in 
Rcf. ^. In the quantum case it is more appropriate to 
rescale the temperature T = T/7 leading to D^ 1 ' — 
rynrpn ^ F(i-fiui, . . . , i-fCj n ). The rescaled bosonic 

Matsubara frequencies u> become continuous variables in 
the extreme quantum limit T — > and because F has a 
small frequency expansion we conclude that the biggest 
contributions to the n phonon loop diagram are of or- 
der 7" as argued in Ref. |^. The diagrammatic expan- 
sion shown in Fig. [l] includes all terms with the mini- 
mal number of internal phonon loops, i.e. includes the 
leading nontrivial order for the electronic self energy and 
the vertex. Note that the leading order in 7 contains 
both, the direct (first term in Fig. |] (b)) and the ex- 
change diagrams fc*-,the vertex. The conventionally de- 
fined ME approacrJl3'til only sums the bubble diagrams 
for the CDW susceptibility, i.e. the exchange diagrams for 
the vertex (which we find to be of .the .same order) are ne- 
glected. Freericks and co-workersli30E3 have considered 
a variety of different vertex corrections to the conven- 
tionally defined ME approximation but they organized 
these in a coupling constant expansion and not in the 
small frequency expansion we considered here. Specific 
vertex corrections^ are of higher order in 7 than the ones 
treated in this work. 

Adding the electron bubbles in Fig. |l| leads to a 
'dressed' phonon propagator X>„ = [(£>") _1 — n n ] _1 with 
a phonon self energy LI n given by 



n„ - (2s + l)T 



E 



Q 'n-\-mQ w 



(14) 



A 



(b) 



> r 

/ n m\ 












FIG. 1. Diagrammatic expansion of the one-particle 

self-energy (a) and of the irreducible vertex of the charge-order 
susceptibility (b) including all terms of order 7. Heavy lines: full 
local Green's function Q n . Dashed lines: free phonon propagator. 
All other possible diagrams are smaller by additional factors of 7 
and are not shown. 



The occurrence of the free phonon propagator restricts 
the bosonic Matsubara frequencies uj n to values on the 
scale of fl. Following the arguments given in Ref. ^ we 
assume that the dominant contributions to H n come from 
frequencies of order t and consequently set the index n 
under the sum to zero.B To leading order in 7 we can 
replace the full by the free Green's function apd. find 
n n = -1/A C + 0{T 2 ) with A c = 3/[4(2s + l)po]M The 
'dressed' phonon propagator reads 



-A 



n 2 



n 2 



(15) 



with the renormalized phonon stiffness and phonon fre- 
quency 



A 



A/(l - A/A c ) 

^(l-A/Ae) 1 / 2 



(16) 
(17) 



Due to this influence of the electrons on the phonons 
the given expansion is actually an expansion in A = Apo 
and 7 = Opo- Maybe a more familiar way to express 
the dressed electron phonon coupling is A = poAT>(0). 
The importance of not neglecting the renormalization of 
phonons by electrons waSj-stressed by seyecal authors, 
both in studies of d = 2t3'E3 and d — cx£J'E3 spatial di- 
mensions. In the limit A —> A c the renormalized phonon 
stiffness A — > 00 and the renormalized phonon frequency 
signaling that the lattice becomes unstable against 
local distortions. In the absence of long range order the 
ground state expectation value of the lattice distortion 
(r) would change from (r) = for A < A c to (r) 7^ for 
A > A c . This polaronic instability has been extensively 
discussed in Ref. ^2] and it is important to note that the 
coupling constant is strongly enhanced and the adiabatic 
parameter is strongly suppressed near the instability. 
Within the small 7 limit we obtain 



F L 

1 nm. — 



Po 

+0 I max 



-(2s + l)A + A 



1 



l + OJ 2 

1 nm 



(18) 



[T 2 A 2 ,7 2 ]) 
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and 



S»i — / Gn 

PO ~f 



(19) 



+0 (max[T 2 /t 2 ,7 2 ]) 



Here, we have rescaled the bosonic frequencies uj nm = 
w nm /0. In Eq. ( JTg| ) as well as in Xn, we can replace 
the full Green's function by the bare one Gn — ► Gn = 
[iu> n — i sgn(u) n )fJLu n + 4t 2 ]/2t 2 to leading order in 7. 
Previous studiesEj asserted that for small phonon fre- 
quencies vertex corrections such as the second term in 
Eq. (|l8| ) can be neglected. We see that this is not quite 
correct: rather, the vertex corrections for 7-^0 affect 



only the diagonal elements of the matrix in Eq. (12), but 
in this diagonal region they are of the same order as the 
contribution from the first part of Eq. (|l8|). 

The importance of the vertex .-corrections has been 
demonstrated in the classical limira. For 7 — > (f2 — > 0) 
we can expand 



1 



1 



l+UJ, 



7,2 

urn 



( 1 ^nm ) * ^n 



(20) 



i.e. the phonon propagator becomes diagonal for 7 = 0. 
In this limit both vertex and self energy corrections to 
the BCS case modify the upper cutoff of the logarith- 
mic divergence in the CDW susceptibility leading to a 
suppression of the CDW transition temperature, which 
has been interpreted rin terms of inelastic scattering of 
electrons by phononsll 

Returning now to the quantum problem we insert 
Eq. (|l8|) into the eigenvalue equation (Il2|) yielding 



(21) 



Po 



-x 1/2 



XT 
Po 



Y - vV 2 e 

^l + Cj 2>Cm 



For generic values of Q/T co (0) we solve Eq. @ numeri- 
cally by restricting the sums to a finite range m < N and 
diagonalizing the resulting matrix exactly numerically. 
Convergence is improved by using an analytical approxi- 
mation for Eq. ( pl| ) for large frequencies uj n > lu^. In the 
large frequency limit the sums in Eq. ( |l9|) and Eq. (|2l]) in- 
cluding the phonon propagator equal their n = m values 
times the factor T^T m 1/(1 -f^) = (fi/2) coth(fi/2T). 
Therefore, the coupling constant (2s+l)A can be replaced 
by the effective coupling constant 



A* = (2s + l)A/[l-(2s + l)A7)] 



with 



Po 

P \n\>N 



( Xn )^ + Cxn/ PQ ) coth(^) 



(22) 
(23) 



where we have introduced x n — ~G n /iu n . For n,m 
0, . . . , N the leading eigenvalue of the matrix 




nAT co (0) 

FIG. 2. CDW transition temperature as function of the typical 
phonon frequency f2 ~ Af — 1 / 2 for spinless electrons. In the main 
panel T co and Q are scaled to the zero frequency value of the CDW 
transition temperature T co (0). Circles connected with dashed lines 
represent results for A = 0.129 but spin s = 1/2. Inset: 
function of £1 on the scale t. 



_ 2X*T 1/2 1/2 AT 1/2 1/2 

J-VJ-nm — An Am An Am 

PO PO 

1 1 



(24) 



1 + [27r(n- m)T/n} 2 1 + [2?r(n + m)T/n} 2 



has to be found. T co is given by the condition that the 
leading eigenvalue equals one. We have checked the con- 
vergence with N; typically we used values ranging from 
N = 256 to N = 4096 yielding a relative error 8T co /T co < 
0.01. For later use we give the weak coupling A — * re- 
sult for the effective coupling constant. In this limit the 
charge density wave transition temperature becomes very 
small and the Matsubara sum can be evaluated as an in- 
tegral yielding 1/A* « l/[(2s + 1)A] - ln(2t/uj N ). 



IV. RESULTS 

The CDW transition tem per ature resulting from a nu- 
merical diagonalization of ( p4| ) is shown in Fig. ^. De- 
creasing the ionic mass, i.e. increasing the phonon fre- 
quency f2, is seen to suppress T co considerably for both 
spinless and spin- 1/2 electrons. The suppression becomes 
very pronounced when Q ~ (2s + l)7rT co (f2 = 0). This is 
in qualitative agreement with the strong suppression of 
the zero temperature order parameter reported in d = 1 
where the CDW state is destroyed when the phonon fre- 
quency is p£ the order of the adiabatic (7 = 0) mean-field 
gap valueo We note that the scaled T co curves fan out at 
larger phonon frequencies: in the extreme quantum limit 
the suppression of T co (f2) relative to T co (Q — 0) depends 
strongly on coupling. This is an indication that quan- 
tum non-Migdal effects (high energy electron self energy 
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£„ w 1/2 Ail 0°) are different from classical non-Migdal 
effects (£„ « AT^°). 

T co (0) increases considerably for spin s = 1/2, essen- 
tially because the spin sum increases the bare suscepti- 
bility by a factor of two. Nevertheless, the change in spin 
s cannot be accounted for by a change of A (i.e. po) only. 
Reducing A by a factor of two when changing the spin 
from s — to s = 1/2 does not bring the transition tem- 
perature curves close to each other as shown in the inset 
of Fig. U (see the circles connected with a solid, s = 0, or 
dashed, s — 1/2, line). The rather different suppression 
of Tcojfoc-s = and s = 1/2 was also found in studies in 
d = l£3'EJ Even on the scale T co (0) the large frequency 
tails are qualitatively different in the two cases. We pro- 
vide some analytic insight in Section V. 

Within the spin-1/2 Holstein model sizable transition 
temperatures might be obtained even in the antiadia- 
batic limit and for modest coupling strengths making the 
CDW transition accessible to QMC simulations. Freer- 
icks and Jarell report a suppression of T co /T co (0)jzt 0.25 
for n/T co (Q) w 7.1 (corresponding to 7 = 0.282).td This 
result, obtained numerically in the antiadiabatic limit, 
indicates an even stronger suppression of T co than ob- 
tained here in the adiabatic limit (7 values used in our 
calculation vary from to 0.08) but is of the same order. 
The comparison of absolute numbers is difficult because 
the work Ref. [l4] does not use the semicircular density of 
states we have chosen. 

It is important to note that T co (0) is not the 
mean-field CDW transition temperature T^ F = 
(8t"f/eTr) exp[— l/(2s + 1)A].lI Even at zero phonon fre- 
quency the charge density wave transition temperature 
is suppressed by thermal fluctuations of the phonon 
fields, e.g. T co (0)/T c ^ F = 0.33,0.40,0.49 for s = and 
A = 0.258, 0.204, 0.156 respectively. In Fig. g T co is plot- 
ted as function of the coupling constant A. The suppres- 
sion of the charge density wave transition temperature 
relative to its 'BCS' or static mean field value due to 
thermal fluctuations (white circle) and due to additional 
quantum fluctuations (white squares) is clearly visible. 
At weak coupling quantum fluctuations of the phonon 
fields are so effective in reducing T co that it seems un- 
likely that CDW materials can be found with a Debye 
temperature 0d much larger than the mean field transi- 
tion temperature T^ F . Higher order corrections in 7 are 
accompanied by higher powers of A and can no longer be 
neglected when the coupling constant is large. The com- 
parison with the full numerical solution of the classical 
DMFT equation given in Ref. || reveals that the expan- 
sion suggested here breaks down already at the surpris- 
ingly small value of 7 A « 10~ 2 . The inset of Fig. || will 
be discussed in Section VI. 

The dependence of the CDW transition temperature 
on the ionic mass is usually discussed in terms of the iso- 
tope exponent a = \ 5P- which is shown in Fig. g. 
Because of the rapidly increasing factor f2/T co the nega- 
tive isotope exponent —a monotonically increases to very 




0.1 0.2 0.3 0.4 0.5 0.6 



X 

FIG. 3. Charge density wave transition temperature as func- 
tion of the coupling strength A for s = 0, calculated from Eq. ( j2l[ ) 
for 7 = (white circles, identical to the modified BCS approxima- 
tion of Ref. 9) and 7 =Ij.05/7r (white squares). Also shown are 
the 'BCS' approximation^ (dashed line) and a numerical solution 
of the classical DMFT equations]^ For A = 0.35 (A7 = 7.6 X 10~ 3 ) 
differences (amplified because T co depends exponentially on param- 
eters) between the exact classical treatment and our approximate 
results become noticeable, suggesting that the expansion breaks 
down. The inset presents data for s = 1/2 and 7 = 0.1/7T. The 
solution of Eq. (Ell) (white squares) is compared to the renor- 
malized Migdal-Eliashberg theory (RME) and to the conserving 
Hartree-Fock approximation (CHF) discussed in Ref. 14. The 7 
expansion can be trusted up to A = 0.18 which corresponds to 
A7 = 7.9 X 10~ 3 . For details sec text. 
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FIG. 4. Isotope exponent as function of Q/t for spinless elec- 
trons. With increasing phonon frequency a drops to huge negative 
values. The circles connected with a thick solid line represents the 
change in a with increasing t for fixed values of A = 0.49 and 
U = 0.015. (ti = 1.65, i 2 = 1.31, ti = 1) The dashed and dotted 
line are extrapolations given by Eq. Q and Eq. (^|), respectively. 
In the latter we have used a cutoff parameter r] = 1.16. The isotope 
exponent for spin-1/2 electrons (dot-dashed line) remains finite. 
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large values for spinless electrons. In the case s = 1/2, 
however, the isotope exponent flattens off as function of 
Sl/t due to the different large frequency behavior of T co . 



V. ANALYTIC EXPRESSIONS 

In the adiabatic and extreme quantum limits we can 
derive explicit analytic expressions for T co . First, we con- 
sider the limit il/T co (Q = 0) — * 0, i.e. we expand (|2l| ) 
around the classical solution obtained in Ref. ^. An ap- 
propriate Ansatz for the eigenvector e n is 



C 



1/2 

Xn 



[l + (\T/ P0 ) Xn \ 



(25) 



Note that, e„ is not quite the classical solution because it 
still depends on 7 via \n- Therefore, after inserting ( p5| ) 
into ( |2l| ) we still have to expand in 7 which finally leads 
to 



1 



2s + 1 



AT v ^ 

77 



Xn 



lass 



(AT/p )x c „ 1 



class 



XK 



r co (o) 



(26) 



with Xn = (V( 2i / W «) 2 + 1 - l)/2t 2 . A" is a numeri- 
cal coefficient independent of temperature but weakly 
dependent on A. Its magnitude is relatively small, e.g. 
K = 1.775 x 1CT 2 for A = 0.156 (s = 0). Q 
The first part of Eq. (|2^) is the classical T co equation.cl 
The sum has a logarithmic T dependence. Therefore, 
expanding T = T co (0) - ST leads to a term [A/T co (0)] ST. 
Consequently, the CDW transition temperature in the 
limit fi/T co (0) < 1 is given by 



T C o(0) 



= 1- K 



T C o(0) 



(27) 



The quadratic drop is clearly seen in Fig. |[ We can use 
Eq. ( p7j ) to derive the isotope exponent in the classical 
limit. We obtain 



n 



Tco(0) 



Tco(0) 



1 



(28) 



The isotope exponent as function of t according to this 
result is depicted in Fig. | as a dashed line. The expo- 
nential decrease of T co (0) with increasing kinetic energy 
t of the electrons is responsible for the sharp upturn of 
the negative isotope exponent. 

In the extreme quantum limit the temperature is not 
a relevant high energy scale. Therefore, the upper cutoff 
frequency introduced above should be given by ujn = T)Q 
to yield results to logarithmic precision with rj being a 
numerical factor of order one. We can replace the phonon 
propagator by 



1 



1 — ~ 2 



(29) 



Moreover, to the same precision all Matsubara sums can 
be replaced by integrals and the temperature occurs only 
as a lower cutoff in the theory. The T co equation ( |2l| ) 
reads simply 



1 = (A* - A) t / ' xH Auj ■ 
Jt 



(30) 



For small frequencies the self energy ( J19|) is given by £„ = 
— iui n A. Therefore, x(uj) = xo( w )/(l + A) which, inserted 
in Eq. (|30|), gives the CDW transition temperature in the 
limit QJT co 3> 1 



T co = r)Q, exp 



1 + A 

A* - A 



,,<> r. M > ( ] ■ CM I 



Note that, f2 enters the calculation as the upper cutoff of 
the integral ( |30| ) and as lower cutoff in the effective inter- 
action strength A* (sec (^2|)) in close resemblance to the 
well known McMillan formula for the superconducting 
transition temper ature.BI Because the effective coupling 
constant A e s depends logarithmically on f2 (via A*) we 
obtain for the isotope exponent in the extreme quantum 
limit 



1 

a= 2 



1 



1 + A 



(1-A/A*) 2 



(32) 



The extreme quantum limit for fixed values of the phonon 
frequency fl is realized for small values of A (or large val- 
ues of t) due to the fast drop of T co (0) with decreasing 
coupling strength (see Fig. f3|). For spinless electrons the 
isotope exponent becomes a w —1/(2 A 2 J 2 ), i.e. a di- 
verges in the limit A — * with power law behavior. The 
divergence of the negative isotope exponent with increas- 
ing kinetic energy of the electrons is shown in Fig. |as a 
dotted line. The isotope exponent for spin-1/2 electrons, 
however, remains finite, equalling a = —3/2 for A — * 0. 
The difference between the cases s = and s = 1/2 is 
due to a near-cancellation of the vertex parts for s = 0. 
Consequently, the leading order of A e ff is 0(X 2 ) for s = 
and 0(A) for s = 1/2. 

In the ME theory for the superconducting transition 
temperature the effect of Coulomb interaction is consid- 
ered by introducing a pseudopotential //*.□ Therefore, it 
is tempting to absorb Coulomb effects on the charge den- 
sity wave transition temperature in a similar way in A*. 
Indeed, a repulsive Coulomb interaction is an electron- 
electron interaction similar to the attractive vertex part, 
first term in Fig. |l|(b), and will change the value of A* sim- 
ply by replacing (2s + 1)A — > (2s + 1)A — U (see also Ref. 
|26|). A* acts as a pseudopotential with a non vanishing 
value for U — 0. Nevertheless, Coulomb interaction will 
also change the electron self energy and the 'vertex cor- 
rections' (second and following terms in Fig. |l|(b)) which 
is a more subtle effect due to the different frequency de- 
pendence of electron-phonon and Coulomb interactions 
and is left for future work. 
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VI. RELATION TO TRADITIONAL 
APPROACHES 

It is traditional to write the CDW transition temper- 
ature in the form T co = W e g exp(— 1/A c ff) defining an 
effective bandwidth and coupling constant. The effective 
parameters W e g and A e ff are usually obtained by fitting 
to experimental data and therefore the functional depen- 
dence on the microscopic properties like electron-phonon 
coupling and adiabatic parameter remains hidden. The 
isotope exponent, however, explores this dependence and 
tests our physical understanding of the CjD.W transition. 

A 'BCS' or static mean field approacbcJ would result 
in T™ F = 2t exp[-l/(2s+l)A], i.e. T co depends exponen- 
tially on A but is independent of 7. For strong coupling 
the 'BCS' theory is clearly wrong because the Holstein 
model in this4imit can be mapped on an effective pseu- 
dospin modeO and T co has to decrease with increasing A. 
Here we have demonstrated that also the weak coupling 
behavior is not predicted correctly and that instead of an 
'BCS' isotope exponent a = an enhanced (s = 1/2) or 
even divergent (s = 0) negative isotope effect is obtained. 

Corrections to the BCS approximation have been con- 
sidered. In the conventional Migdal-Eliashberg (ME) 
theorya of electron-phonon effects in solids electron self 
energies arising from the interaction with the phonons are 
included but all vertex corrections are neglected. Other 
authorstSEJ have argued that one should include also the 
feedback of the electrons on the phonon system but still 
neglecting the vertex corrections, thus defining a renor- 
malized Migdal-Eliashberg theory (RME). In our nota- 
tion this corresponds to dropping the second (A) term 
in Eq. (|2l]). The RME approximation predicts an in- 
correct isotope effect because T co depends on ion mass 
only via the dependence of Xm on the normal state self 
energy, and comparison to our results shows that this de- 
pendence is not the leading one. At very weak coupling 
the RME theory reduces to the BCS theory (up to log- 
arithmic corrections) and predicts a — 0. Freericks and 
co-woijkers added extra vertex corrections to the RME 
theoryEj but these are higher order in 7 then the effects 
we studied in this paper. 

We have plotted the charge density wave transition 
temperature for various coupling strengths A and fixed 
adiabatic parameter 7 = 0.1 /ir and s = 1/2 in the inset 
of Fig. H for both the RME approach and the present 
work. The charge density wave transition temperature 
obtained within RME theory is always larger than the 
one obtained from the systematic 7 expansion presented 
here. A similar plot for a larger value of 7 = 0.5 / y/n 
has been presented in Ref. Q where the RME results 
were compared with QMC data, showing that the RME 
theory indeed overestimates the transition temperature 
for weak coupling strength at least in the antiadiabatic 
limit. Note that, the RME results presented in Ref. |lj 
are not precisely identical to ours because we have in- 
cluded the phonon self energy only to leading order in 7, 



but we believe that this does not effect the conclusions 
drawn here. It is well known that T co (A) is maximized 
at a A max < 00. [T co — > as A — * because the cou- 
pling vanishes; T co — > as A — ► 00 because in this limit 
a picture of almost completely localized particles with 
intersite interactions ~ t 2 /X applies.] The present ap- 
proach underestimates both A max and T" max , as may be 
seen (Fig. ||) by comparison to the classical treatment 
of Ref. H (to which the present treatment reduces as T co 
becomes greater than fi) and by comparison to QMC 
data (Ref. |l4|) and indeed is not expected to apply for 
A 00 A max . 

A theory neglecting the renormalization of the phonon 
propagator by electrons predicts that T co is maximized 
at A = 00. As an example we have plotted in the in- 
set of Fig. |^ the conserving .Hartree-Fock approximation 
discussed by Freericks et a/.Ej This corresponds to using 
Eq. ( pl| ) but with A replaced by A. In the weak coupling 
limit the difference between A and A is negligible, so the 
CHF theory gives a good estimate of the CDW transition 
temperatures. The comparison with QMC data given in 
Ref. 14 demonstrate that the CHF theory, and there- 
fore the theory given here, can (at least for s = 1/2 and 
weak coupling strengths) be extended to the antiadia- 
batic limit. 

Finally, we point out that the systematic expansion 
in 7 discussed in this work reproduces the results of 
the renormalized ME theory when applied to the cal- 
culation of the superconducting transition temperature 
T c . Following Freericks and ScalapinoE3 we find T c w 
r/flexp[—(l + A)/A] in the extreme quantum limit. For 
the half-filled Holstein model the charge density wave 
transition temperature is always larger than the super- 
conducting transition temperature if Q < 00; only as 
ft — > 00 do they become equal.E3 The superconducting 
transition temperature is controlled by a logarithmic di- 
vergence with upper cutoff set by the phonon frequency 
whereas the corresponding upper cutoff in the theory 
for T co is enhanced towards 2i by the first diagram in 
Fig. 0(b). 



VII. CONCLUSION AND COMPARISON TO 
EXPERIMENT 

We have used the dynamical mean field method to ob- 
tain a systematic and tractable treatment of the isotope 
effect on charge density wave transitions. Our method is 
applicable to the physically relevant case in which both 
phonon frequency and transition temperature are small 
relative to the electronic bandwidth, and in combination 
with the results of Ref. || may be generalized to the cru- 
cial case in which other interactions beyond the electron- 
phonon one are important. 

Quantum phonon effects are found to suppress the 
charge density wave transition temperature severely. 
This is in agreement with previous work on one dimen- 
sional models demonstrating that the suppression of T co 
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is not due to the strong quantum fluctuations charac- 
teristic of one spatial dimension. The isotope exponent 
of T co for spin- 1/2 electrons is enhanced compared to 
the ME prediction. In the case of spinless electrons the 
isotope exponent may even be divergent. Although be- 
yond the scope of the theory, the consideration of phonon 
self energy effects lead to an optimal charge density wave 
transition temperature for values of the coupling strength 
slightly below the polaronic instability at A c but peak po- 
sition and height are underestimated. 

Our results suggest that it will be difficult to find CDW 
systems with mean field transition temperatures much 
less than typical Debye frequencies. It is of interest to ex- 
amine data with this in mind. For quasi-one-dimcnsional 
CDW materials the mean field CDW transition tem- 
perature is typically larger than room temperature and 
therefore larger than or of the same order as the De- 
bye temperature 9d- In TaS3, K0.3M0O3 and (TaSe^I 
one deduces experimentally values of T C ^ F p. 540 A, 
T™ F w 330 and T™ F « 1040A respectively^ Quasi- 
two-dimensional CDW materials include the dichalco- 
genides and the molybdenum bronzes and oxides. For 
these materials it is more difficult to deduce mean field 
values of the CDW transition temperature because the 
CDW transition is not accompanied by a metal insula- 
tor transition. For 2H — TaSe2 one might argue from 
measurements of the optical conductivityEa for a value of 
T^ F 580A' similar to the quasi-one-dimensional mate- 
rials. Charge order transitions are also observed in many 
two and three dimensional transition metal oxides includ- 
ing some members of the high-temperature superconduc- 
tor and 'CMR' family.Ej Frequently, the Coulomb repul- 
sion of the charge carriers is responsible for the charge or- 
der transition and therefore beyond the scope of this pa- 
per. Nevertheless, the phenomenology is very similar and 
e.g. for Lai.67Sr .33NiO4 one might deduce T^ F = 750A. 
As was suggested by our theory the case of Od 3> T^ F 
seems to be extremely rare. 

Our results bear at least a qualitative similarity to 
data obtained on 'CMR' manganites. Recently, a colos- 
sal isotope shift of the CDW transition temperature in 
Ndo.5Sr . 5 Mn03 was reportedJa With magnetic field the 
negative isotope exponent varies from ~ 1.2 to ~ 4.5. 
Within the double exchange picture of the manganites 
an increase of the applied magnetic field means an in- 
crease of the kinetic energy t. Due to the large Hunds 
rule coupling of the conduction electron spins to the core 
spins of the Mn ions the conduction electrons are fully 
spin polarized. Therefore we should compare to the case 
s = 0. To demonstrate the effect of an increasing t in our 
calculation we draw a line for fixed A and Q in Fig. ||. 
Indeed, we find a monotonically increasing negative iso- 
tope effect with increasing t. 
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